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Abstract 

Background: One-carbon metabolism is the basement of nucleotide synthesis and the methylation of DNA linked to cancer 
risk. Variations in one-carbon metabolism genes are reported to affect the risk of many cancers, including renal cancer, but 
little knowledge about this mechanism is known in Chinese population. 

Methods: Each subject donated 5 mL venous blood after signing the agreement. The study was approved by the 
Institutional Review Board of the Nanjing Medical University, Nanjing, China. 18 SNPs in six one-carbon metabolism-related 
genes (CBS, MTHFR, MTR, MTRR, SHMT1, and TYMS) were genotyped in 859 clear cell renal cell carcinoma (ccRCC) patients 
and 1005 cancer-free controls by the Snapshot. 

Results: Strong associations with ccRCC risk were observed for rs706209 (P = 0.006) in CBS and rs9332 (P = 0.027) in MTRR. 
Compared with those carrying none variant allele, individuals carrying one or more variant alleles in these two genes had a 
statistically significantly decreased risk of ccRCC [P = 0.001, adjusted odds ratio (OR) = 0.73, 95% confidence interval 
(CI) = 0.06-0.90]. In addition, patients carrying one or more variant alleles were more likely to develop localized stage disease 
(P = 0.002, adjusted OR =1.37, 95%CI = 1.11-1 .69) and well-differentiated ccRCC (P<0.001 , adjusted OR = 1 .42, 95%CI = 0.87- 
1.68). In the subgroup analysis, individuals carrying none variant allele in older group (P = 0.007, adjusted OR = 0.67, 
95%CI = 0.49-0.91), male group (P = 0.007, adjusted OR = 0.71, 95%CI = 0.55-0.92), never smoking group (P = 0.002, adjusted 
OR = 0.68, 95%CI = 0.53-0.88) and never drinking group (P<0.001, adjusted OR = 0.68, 95%CI = 0.53-0.88) had an increased 
ccRCC risk. 

Conclusions: Our results suggest that the polymorphisms of the one-carbon metabolism-related genes are associated with 
ccRCC risk in Chinese population. Future population-based prospective studies are required to confirm the results. 

Citation: Zhang L, Meng X, Ju X, Cai H, Li P, et al. (2013) One-Carbon Metabolism Pathway Gene Variants and Risk of Clear Cell Renal Cell Carcinoma in a Chinese 
Population. PLoS ONE 8(11): e81129. doi:10.1371/joumal.pone.0081 129 

Editor: Yong-Gang Yao, Kunming Institute of Zoology, Chinese Academy of Sciences, China 

Received April 26, 2013; Accepted October 9, 2013; Published November 21, 2013 

Copyright: © 2013 Zhang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was supported by the Program for Development of Innovative Research Team in the First Affiliated Hospital of Nanjing Medical University, 
Provincial Initiative Program for Excellency Disciplines of Jiangsu Province, the National Natural Science Foundation of China [grant number 81 1 71 963, 81 1 02089, 
and 81201998]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 

* E-mail: qinchao@njmu.edu.cn (CQ); drcjyin@gmail.com (CY) 

9 These authors contributed equally to this work. 



Introduction 

Renal cell carcinoma (RCC), one lethal urologic cancer 
accounting for 2% of all cancer deaths, is the third most common 
cancer of the genitourinary tract, just next to prostate and bladder 
cancer [1]. The most common subtype of RCC is the clear cell 
type (ccRCC), accounting for 75% of all cases. The incidence of 
RCC has increased rapidly worldwide over the past few decades 
[2] . Although the exact causes of RCC remain unknown, a few 
well-established lifestyle risk factors have been identified, includ- 
ing: cigarette smoking, obesity, hypertension and diabetes [3]. 
Recent studies have demonstrated that genetic polymorphisms of 
candidate genes were associated with susceptibility and/or 
prognosis of cancers, including RCC [4,5,6,7]. 



One-carbon metabolism refers to a system of interdependent 
metabolic pathways facilitating the transfer of one-carbon units, 
which are the basement of nucleotide synthesis and the 
methylation of DNA. Methylation is responsible for gene 
expression control, chromatin structure stability and the mainte- 
nance of genomic stability. It has been proved that tumor- 
suppressor genes are silenced by hypermethylation of the promoter 
region in various cancers, leading to carcinogenesis and cancer 
progression [8]. In RCC, genes encoding von Hippel-Lindan 
(VHL), E-cadherin and RASSF1A have been reported to be 
inactivated by promoter hypermethylation [9,10,11]. DNA 
hypomethylation begets chromosomal instability as a result of 
changes in chromatin structure [12,13]. 
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The connection between the polymorphisms in one-carbon 
metabolism pathway genes, such as methylenetetrahydrofolate 
reductase (MTHFR), has been widely studied in various cancers 
[14,15,16,17,18]. The MTHFR promoter polymorphism 
rs3737965 has been intensively studied, but the results were 
conflicting, especially in lung cancer [16,19,20,21,22]. Besides, two 
large case-control studies for RGC have been conducted in Europe 
to analyze the candidate genes in one-carbon metabolism 
pathway. The studies found the association between the one- 
carbon metabolism pathway genes polymorphisms and RGC risk 
[14,15]. However, there are some discrepancies in these two 
studies, as evidenced by the confusing connection between the 
TYMS polymorphisms and RCC risk. Given the importance of this 
pathway in carcinogenesis, ethnic differences and the discrepan- 
cies in previous studies, the present case-control study was 
conducted to define the association between the one-carbon 
metabolism pathway genes polymorphisms and the ccRCC risk in 
Chinese population, through testing the polymorphism in six key 
genes {CBS, MTHFR, MTR, MTRR, SHMT1, and TYMS) in this 
pathway. 

Materials and Methods 

Study Subjects 

859 ccRCC patients and 1005 cancer-free controls were 
enrolled in our case-control study. All subjects were genetically 
unrelated ethnic Han Chinese recruited from May 2004 to 
January 2012 at the First Affiliated Hospital of Nanjing Medical 
University, Nanjing, China, the Affiliated Hospital of Medical 
College Qingdao University, Qingdao, China and the Third 
Affiliated Hospital of Soochow University, Changzhou, China 
[23]. The cases were newly diagnosed with incident ccRCC and 
had been histopathologically confirmed without history of other 
malignancies and previous chemotherapy or radiotherapy. A 
standard questionnaire was administered through face-to-face 
interviews by trained interviewers to collect demographic data and 
related information. Each subject donated 5 mL venous blood 
after signing the agreement. The disease was classified by the 
World Health Organization criteria and stage according to the 
American Joint Committee on Cancer (AJCC) TNM (tumor- 
node-metastasis) classification. The Fuhrman scale was used to 
assess tumor nuclear grade and the disease stage was divided into 
two subgroups (localized group: stage I and stage II and advanced 
group: stage III and stage IV). The controls were individuals 
without history of cancer who were seeking health care in the 
outpatient departments matched to the cases' sex and age (±5 
years) on frequency. We subdivided the patients and controls into 
two groups (younger group: age ^57 and older group: age >57) 
according to the mean age of the cohort (57 years old). Smokers 
were defined as those who smoked daily for more than 1 year. 
Drinkers were those who drank at least 3 times per week and more 
than 6 months. The study was approved by the Institutional 
Review Board of the Nanjing Medical University, Nanjing, China. 

DNA extraction and genotyping 

Genomic DNA was seperated and purified from leucocytes of 
peripheral blood by proteinase K digestion and phenol/chloro- 
form extraction. Genotyping of the polymorphism was performed 
by the Snapshot. The SNaPshot SNP assay was performed to 
detect the dimorphism at the two SNP loci. PCR was performed in 
a 20 n.1 reaction mixture containing 1 mJ (10 ng) of template DNA, 
1 uM of each primer, 0.3 inM of each deoxynucleotide triphos- 
phate, 3.0 mM of MgC12, and 1 U HotStarTaq polymerase 
(Qiagen Inc., USA) with lx HotStarTaq buffer. The PCR 



program consisted of an initial melting step of 15 minutes at 95°C; 
11 cycles of 20 seconds at 94°C, 40 seconds at 65°C-0.5°C/cycle, 
and 90 seconds at 72°C; 24 cycles of 20 seconds at 94°C, 30 
seconds at 59°C, and 90 seconds at 72°C; and a final elongation 
step of 2 minutes at 72°C. To purify the PCR products, 1 U SAP 
and 1 U Exonuclease I were mixed with 10 PCR product for 
1 hour at 37°C and 15 minutes at 75°C.The extension reaction 
was performed in a 10 (0.1 reaction mixture containing 5 JU.1 of the 
SNaPshot Multiplex Kit (Applied Biosystems, USA), 2 ul of 
purified PGR products, 0.8 |J,M of the extension reaction primer, 
and 2 JU.1 water. The PCR program was 1 minute at 96°C; 28 
cycles of 10 seconds at 96°C, 5 seconds at 50°C, and 30 seconds at 
60°C; and 4°C as the holding temperature. Finally, 10 ul of the 
extension product was purified with 1 U SAP for 1 hour at 37°C 
and inactivated for 15 minutes at 75°C.The resulting data were 
analyzed with an ABI3130XL sequencer and GeneMapper™ 4.0 
Software (Applied Biosystems, Co. Ltd., USA). All the sequences 
of primers for each SNP are listed in Table SI. Positive controls by 
sequencing and negative controls were used to ensure the accuracy 
of genotyping. To ensure quality control (QC), genotyping was 
performed by researchers blinded to the case/control status of the 
subjects, and a random sample of 5% of the cases and controls was 
genotyped twice by different persons, with a reproducibility of 
100%. 

Statistical analysis 

Statistical analysis were performed by software SAS 9.1.3 (SAS 
Institute, Cary, NC). Allele frequencies were tested against 
departure from Hardy- Weinberg equilibrium through goodness- 
of-fit x2 test before analysis. Differences in the distributions of 
demographic characteristics, selected variables, and frequencies of 
genotypes between cases and controls were evaluated by Student's 
t-test (for continuous variables) or %2-test (for categorical 
variables). OR and 95% CIs from unconditional logistic regression 
analysis with the adjustment for possible confounders were used to 
estimate the association between the polymorphism and the risk of 
ccRCC. All of the statistical tests were two-sided and P<0.05 was 
considered statistically significant. 

Results 

Characteristics of ccRCC patients and controls 

Clinical and pathologic characteristics of the 859 cases and 1005 
controls are presented in Table 1 (negative results listed in Table 
S2). No significant differences were found between the cases and 
controls in regards to age, gender and drinking status (all i*>0.05). 
However, there were more high BMI, smoking habit, hypertension 
and diabetes subjects in ccRCC patients than those in controls 
(P= 0.003, P= 0.040, P<0.001 and P= 0.040, respectively). 
Among 859 ccRCC cases, 577 (67.2%) patients were diagnosed 
with stage I disease, 170 (19.8%) patients with stage II, 57 (6.7%) 
patients with stage III, and 55 (6.4%) patients with stage IV. The 
percent of nuclear grade from I to IV was 20.5%, 54.7%, 19.9%, 
and 4.9%, respectively. 

Association of SNPs with the risk of ccRCC 

As shown in Table 2, except the rs 1788484 in CBS gene and 
rsl0520873 in MTRR gene, the genotype distributions of all other 
SNPs were in agreement with Hardy- Weinberg equilibrium in 
controls. Therefore, the rs 1788484 and rs 1052087 3 were excluded 
for further analysis. Table 2 shows ccRCC risks were associated 
with each single polymorphism. We presented each of the 
associations in additive, dominative and recessive models. In 
CBS gene, the genotype of rs706209 polymorphism distribution 
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Table 1. Distribution of selected variables between the clear 
cell renal cell carcinoma cases and control subjects. 





Variables 


Cases (n 


= 859) 


Controls 
(n = 1005) 


lvalue" 




N 


% 


N % 




Age (years) (mean ± SD)57.0±11.7 




57.2±12.4 


0.758 


BMI (kg/m2) (mean 
SD) 


± 24.2 ±2.9 




23.8±3.2 


0.003 


Gender 


Male 


557 


64.8 


688 68.5 


0.098 


Female 


302 


35.2 


317 31.5 




Smoking status 


Never 


530 


61.7 


666 66.3 


0.040 


Ever 


329 


38.3 


339 33.7 




Drinking status 


Never 


619 


72.1 


757 75.3 


0.110 


Ever 


240 


27.9 


248 24.7 




Hypertension 


No 


522 


60.8 


675 67.2 


<0.001 


Yes 


337 


39.2 


251 25.0 




Diabetes 


No 


748 


87.6 


870 86.6 


<0.001 


Yes 


111 


12.4 


56 5.6 




Clinical stage 


I 


577 


67.2 






II 


170 


19.8 






III 


57 


6.7 






IV 


55 


6.4 






Grade 


I 


176 


20.5 






II 


470 


54.7 






III 


171 


19.9 






IV 


42 


4.9 






a Student's t-test for age and BMI distributions between cases and controls; two- 



sided x2-test for others selected variables between cases and controls. 
doi:1 0.1 371 /journal.pone.0081 1 29.t001 

between cases and controls was significantly different (P= 0.006). 
In the dominant model, when compared with the homozygous 
wild-type reference group, TC+CC genotype was significantly 
associated with a decreased risk of ccRCC (P= 0.002, adjusted 
OR=0.76, 95% CI = 0.62-0.91). InMTXflgene, the genotype of 
rs9332 polymorphism distribution between cases and controls was 
significantly different (P= 0.027). In the dominant model, CT+TT 
genotype was significantly associated with a decreased risk of 
ccRCC {P= 0.010, adjusted OR = 0.74, 95%CI = 0.60-0.92), 
when compared with the CC group. The rs2966952 seemed to 
be a bordering positive result. In the recessive models, TT 
genotype was slighdy associated with an increased risk of ccRCC 
(P=0.050, adjusted OR= 1.34, 95%CI = 1.00-1.80), compared 
with the CC+CT genotype. In TTMS gene, the genotype of 
rs2853741 polymorphism distribution between cases and controls 
had no significantly difference (P= 0.096). In the dominant model, 
TC+CC genotype was significantly associated with an increased 
risk of ccRCC (P= 0.032, adjusted OR = 1.25, 95%CI=1.02- 
1.53), when compared with the TT group. Similarly, for the 



rs699517, the polymorphism distribution between cases and 
controls had no significandy difference (P= 0.108). In the recessive 
models, CC genotype was seemed to be associated with an 
increased risk of ccRCC (P=0.172, adjusted OR =1.25, 
95%CI= 1.03-1.51). 

Combined Analysis of Polymorphisms 

Rs706209 or rs9332 alone being associated with ccRCC risk, 
we combined these two polymorphisms based on the number of 
the variant alleles. As shown in Table 3, statistically significant 
difference was obviously found between the subgroups and the risk 
of ccRCC. Meanwhile, individuals carrying 1, 2, 3, or 4 variant 
alleles in these two genes were associated with a statistically 
significandy decreased risk of ccRCC compared with individuals 
carrying none variant alleles (P— 0.001, adjusted OR = 0.73, 
95%CI = 0.06-0.90). 

Stratified Analysis of the Two Polymorphisms and 
Clinicopathologic Characteristics and Risk of ccRCC 

We further evaluated the association between the combined 
genotypes of rs706209 and rs9332 polymorphisms and clinico- 
pathologic characteristics of ccRCC. These two polymorphisms 
were combined based on the number of the variant alleles. As 
shown in Table 4, a significantly increased risk appeared in 
ccRCC patients with localized stage (P= 0.002, adjusted 
OR= 1.37, 95%CI= 1.11-1.69) and in patients with well-differ- 
entiated ccRCC (P<0.001, adjusted OR= 1.42, 95%CI= 1.14- 
1.76). 

As shown in Table 5, the analysis was stratified by age, gender, 
smoking status, and drinking status. Individuals carrying none 
variant allele in older group (P= 0.007, adjusted OR = 0.67, 
95%CI = 0.49-0.91), male group (P=0.007, adjusted OR = 0.71, 
95%CI = 0.55-0.92), never smoking group (P= 0.002, adjusted 
OR =0.68, 95%CI = 0.53-0.88) and never drinking group 
(P<0.001, adjusted OR = 0.68, 95%CI = 0.53-0.88) had an 
increased ccRCC risk. 

Discussion 

Our study supports a role that CBS, MTRR, and TTMS play in 
modifying ccRCC risk. Meanwhile, we also notice the discrepan- 
cies between our findings and the results reported by the two 
European studies. Moore et al. reported that MTHFR and TTMS 
had effect on RCC risk, while CBS, MTR, or MTRR had no [14]. 
Gibson et al. reported that the strongest association between RCC 
risk and MTHFR, not CBS, MTR, MTRR, SHMT1 or TTMS, was 
observed [15]. The discrepancies might be due to the different 
aspect of study, different ethnic population and other unknown 
factors. Our current study is focus on the ccRCC, which is one 
part of RCC. The pathologic characteristics of ccRCC in our study 
may be relatively simple, when compared with these two 
European articals. In the dbSNP database, the minor allele 
frequency (MAE) of rs706209 was T allele = 0.403 in European, 
but the MAF in Chinese was C allele = 0.402, close to our data C 
allele = 0.482. As involved in one-carbon metabolism pathway, 
the CBS gene encodes cystathionine beta synthase (CBS), which is 
the central enzyme in the transsulfuration pathway irreversibly 
metabolizing homocysteine (Hey) to cystathionine [24] . Evidences 
showed that functional CBS gene SNP could impair CBS gene 
function, leading to an increase in the concentration on tHcy that 
further influenced aberrant DNA methylation patterns and CBS 
gene SNP had associations with lung cancer, colorectal cancer and 
head and neck squamous cell carcinoma [16,25,26]. Whether 
rs706209 has the similar influence on the CBS gene function needs 
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Table 2. Polymorphism 


in one-carbon metabolism 


genes and clear 


cell renal eel 


carcinoma risk. 








Position 


MAF Cases, n (%) 


Controls, n (%) 


/'-value' 


Adjusted OR (95% Cl) b 


P trend* 


HWE 


CBS 


rsl 788484 promoter 


0.222 












CC 


535(62.3) 


543(54.0) 


<0.001 


Ref. 


0.022 


<0.001 


CT 


275(32.0) 


425(42.3) 




0.62(0.51-0.76) 






TT 


49(5.7) 


37(3.7) 




1 .27(0.79-2.20) 






TT vs. CC+CT 






0.038 


1.53(0.97-2.41) 






CT+TT vs. CC 






<0.001 


0.67(0.55-0.82) 






rs706209 3'UTR 


0.329 












TT 


401(46.7) 


396(39.4) 


0.006 


Ref. 


0.007 


0.379 


TC 


358(41.7) 


480(47.8) 




0.76(0.62-0.92) 






CC 


100(11.6) 


129(12.8) 




0.74(0.55-1.01) 






CC vs. TT+TC 






0.434 


0.87(0.65-1.16) 






TC+CC vs. TT 






0.002 


0.76(0.62-0.91) 






MTRR 


rs2966952 promoter 


0.333 












CC 


371(43.2) 


463(46.1) 


0.120 


Ref. 


0.064 


0.978 


CT 


374(43.5) 


438(43.6) 




1.07(0.88-1.31) 






TT 


114(13.3) 


104(10.4) 




1.40(1.03-1.91) 






TT vs. CC+CT 






0.050 


1.34(1.00-1.80) 






CT+TT vs. CC 






0.213 


1.13(0.93-1.36) 






rs9332 3'UTR 


0.188 












CC 


651(75.8) 


708(70.5) 


0.027 


Ref. 


0.007 


0.434 


CT 


1 96(22.8) 


275(27.4) 




0.75(0.60-0.93) 






TT 


12(1.4) 


22(2.2) 




0.65(0.31-1.33) 






TT vs. CC+CT 






0.203 


0.70(0.34-1 .43) 






CT+TT vs. CC 






0.010 


0.74(0.60-0.92) 






rs10520873 3'UTR 


0.311 












AA 


593(69.0) 


648(64.5) 


0.063 


Ref. 


0.098 


0.032 


AG 


240(27.9) 


331(32.9) 




0.79(0.64-0.97) 






GG 


26(3.0) 


26(2.6) 




1.15(0.65-2.04) 






GG vs. AA+AG 






0.566 


1.24(0.70-2.19) 






AG+GG vs. AA 






0.038 


0.82(0.67-1.00) 






TYMS 


rs2853741 promoter 


0.411 












TT 


264(30.7) 


356(35.4) 


0.096 


Ref. 


0.046 


0.982 


TC 


439(51.1) 


484(48.2) 




1.25(1.01-1.54) 






CC 


156(18.2) 


165(16.4) 




1.24(0.93-1.64) 






CC vs. TT+TC 






0.321 


1.08(0.84-1.39) 






TC+CC vs. TT 






0.032 


1.25(1.02-1.53) 






rs699517 3'UTR 


0.3 












TT 


420(48.9) 


537(53.4) 


0.108 


Ref. 






TC 


358(41.7) 


391(38.9) 




1.25(1.02-1.52) 


0.035 


0.619 


CC 


81(9.4) 


77(7.7) 




1.26(0.89-1.78) 






CC vs. TT+TC 






0.172 


1.25(1.03-1.51) 






TC+CC vs. TT 






0.051 


1.15(0.82-1.61) 







a Two-sided x2-test for either genotype distributions or trend between the cases and controls. 

b genotype-specific ORs were adjusted for age, gender, BMI, smoking status, drinking status, diabetes and hypertension in logistic regression model. 
doi:1 0.1 371 /joumal.pone.0081 1 29.t002 
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Table 3. Combined genotype 
carcinoma. 


frequencies a 


mong the case patients and control subjects and their collection to clear cell renal cell 




Variables 


Case (n = 


859) 


Controls (n 


= 1005) 


lvalue 3 


Adjusted OR (95% Cl) b 




n 


% 


n 


% 






rs706209 and rs9332 


Number of risk allele 


0 


304 


35.4 


286 


28.5 






1 


358 


41.7 


432 


43.0 


0.022 


0.79(0.64-0.99) 


2 


171 


19.9 


237 


23.6 


0.003 


0.68(0.52-0.88) 


3 


26 


3.0 


49 


4.9 


0.006 


0.48(0.29-0.81) 


4 


0 


0 


1 


0 






Recombined groups 


0 


304 


35.4 


286 


28.5 






1-4 


555 


64.6 


719 


71.5 


0.001 


0.73(0.06-0.90) 


a Two-sided yT. test for the distributions of genotypes. 

b Adjusted for age, gender, BMI, smoking status, drinking status, hypertension 


and diabetes 


in logistic 


regression model. 





doi:1 0.1 371 /journal.pone.0081 1 29.W03 



further functional studies. MTRR catalyzes reductive methylation 
of cob (Il)alamin by using SAM as a methyl donor to reactivate 
MTR. Thus, MTRR may act as a key regulator of the 
homocysteine conversion to methionine [27]. Our study demon- 
strated that CT+TT genotype of rs9332 was significantly 
associated to a decreased risk of ccRCC. Other report demon- 
strated that rs9332 in MTRR had connections with spina bifida 
and conotruncal heart defects [28]. Although litde information 
about the function of this polymorphism was known, it is known 
that 3'UTR of the MTRR gene may influence miRNA or siRNA 
binding target site, causing the degradation of the mRNA, or 
inhibiting translation initiation. This polymorphism has a protec- 
tive role in the risk of ccRCC. More researches should be done to 
reveal the relationship between them. 

Moreover, we found that TYMS SNPs (rs2966952 and 
rs699517) had some of connections with ccRCC risk. Skibola 
et al. reported that C allele in rs699517 was approximately 100% 
correlated with TYMS 1496 insertion and T allele in rs699517 was 



also approximately 100% correlated with TYMS 1496 deletion. 
The TYMS 1496 deletion polymorphism was associated with 
decreased mRNA stability and low expression in tumor tissue than 
the wild type polymorphism [29,30,31]. Evidence demonstrated 
that the level of thymidylate synthetase (TS) activity was correlated 
with both the progression of the stage and the increase of the grade 
of RCC, and the activity of TS was approximately 5-fold higher in 
RCC compared with normal kidney [32]. 

There were several limitations in our present study. Firstly, our 
standard questionnaire contained no environmental factors, such 
as occupational exposure, and personal habits, such as diet and 
physical activity. So our study may have a limited statistical power, 
such as that we could not do further analysis in gene-environment 
interaction. Secondly, more and more studies indicated that gene- 
gene interaction may also have contribution to the risk of cancer 
[33,34]. In the current study, we did litde analysis about the 
interaction among genes in the one-carbon metabolism pathway. 
So the results may exist some interaction bias. Furthermore, our 



Table 4. Association between the combined genotypes of rs706209 and rs9332 polymorphisms and clinicopathologic 
characteristics of clear cell renal cell carcinoma. 



Variables Risk allele />-value" Adjusted OR (95% Cl) b 





0 


1-4 








n % 


n % 






Control (n = 1005) 


286(28.5) 


719(71.5) 




Ref. 


Case (n = 859) 


Clinical stage 


Localized (l+ll) 


265(35.5) 


482 (64.5) 


0.002 


1.37 (1.11-1.69) 


Advanced (lll+IV) 


39 (34.8) 


73 (65.2) 


0.160 


1.30 (0.86-1.98) 


Grade 


Well- differentiated (l+ll) 


234 (36.2) 


412 (63.8) 


<0.001 


1.42 (1.14-1.76) 


Moderately and Poorly differentiated(lll+IV) 


70 (32.9) 


143 (67.1) 


0.199 


1.21 (0.87-1.68) 



a Two-sided -/2 test for the distributions of genotypes. 

b Adjusted for age, gender, BMI, smoking status, drinking status, hypertension and diabetes in logistic regression model. 
doi:1 0.1 371 /journal.pone.0081 1 29.t004 



PLOS ONE | www.plosone.org 



5 



November 2013 | Volume 8 | Issue 11 | e81 1 29 



One-Carbon Metabolism Genes Relate to ccRCC Risk 



Table 5. Association between the combined genotypes of rs706209 and rs9332 polymorphisms and clear cell renal cell carcinoma 
in stratified analysis. 





Variables 


Risk allele 








P-value" 


Adjusted OR (95% Cl) b 




0 




1-4 










Case(n,%) 


control(n,%) 


case(n,%) 


control(n,%) 






Age 


< = 57 


161(36.2) 


164(30.4) 


284(63.8) 


376(69.6) 


0.054 


0.78(0.60-1 .02) 


>57 


143(34.5) 


122(26.2) 


271(65.5) 


343(73.8) 


0.007 


0.67(0.49-0.91) 


Gender 


Male 


192(34.5) 


188(27.3) 


365(65.5) 


500(72.7) 


0.007 


0.71(0.55-0.92) 


Female 


121(37.1) 


98(30.9) 


190(62.9) 


219(69.1) 


0.105 


0.77(0.55-1 .08) 


Smoke 


Never 


193(36.4) 


187(28.1) 


337(63.6) 


479(71.9) 


0.002 


0.68(0.53-0.88) 


Ever 


111(33.7) 


99(29.2) 


218(66.3) 


240(70.8) 


0.207 


0.83(0.59-1.17) 


Drink 


Never 


229(37.0) 


202(26.7) 


390(63.0) 


555(73.3) 


<0.001 


0.68(0.53-0.88) 


Ever 


75(31.3) 


84(33.9) 


165(68.8) 


164(66.1) 


0.537 


0.83(0.59-1.17) 



a Two-sided yT. test for the distributions of genotypes. 

b Adjusted for age, gender, BMI, smoking status, drinking status, hypertension and diabetes in logistic regression model. 
doi:1 0.1 371 /journal.pone.0081 1 29.t005 



study was designed as a hospital-based study, so the possibility of 
selection bias of subjects could not be ruled out. Thirdly, the 
genotype distributions of rsl788484 in CBS gene and rsl0520873 
in MTRR gene were departed from Hardy- Weinberg equilibrium 
in controls. The probability of genotyping error was very low. 
Meanwhile, Galbiatti et al. also reported that polymorphism of 
CBS gene was not in Hardy- Weinberg equilibrium [26]. The 
departure from the Hardy- Weinberg equilibrium may result from 
selection bias, disease model adopted, and evolutionary factors 
which may influence changes in the genotype frequencies [35,36] . 
On the other hand, this disequilibrium should be expected, in the 
case that it reflected biological and genetic characteristics in 
complex disease models [37]. 

In summary, the data indicates that the common variation in 
CBS, MTRR and TYMS may significantly modify ccRCC risk. 
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